High-field transport properties of bulk Si: A test for the Fokker-Planck approach 
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High electric-field transport parameters are calculated using an analytical Fokker-Planck ap- 
proach (FPA), where transport is modeled as a drift-diffusion process in energy space. We have 
applied the theory to the case of Si, taking into account the six intervalley phonons, aiming to 
test the FPA. The obtained results show a quite reasonable agreement with experimental data and 
Monte Carlo simulations confirming in this case that the FPA works very well for high enough 
electric fields. 
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More than three decades ago the Fokker-Planck approach (FPA) was proposed as an alternative for the Boltzmann 
transport equation (BTE) in the calculation of semiconductor transport properties [jl]Q] . Recently the theory has been 
revisited in a series of works where a relatively detailed analysis of both mathematical and physical aspects of 
, this formalism was developed. In these papers the same model system was considered in calculations using the FPA 
and the BTE by means of the Monte Carlo method. The results of both approaches showed good agreement in the 
high electric-field regime for the mentioned model system. In spite of this, there remains a certain degree of doubt 
about how the FPA could handle a more realistic model of a concrete semiconductor with several possible scattering 
mechanisms and complicated band structure. The FPA considers transport, in the opposite regime of the ballistic 
one, as a certain diffusive-drifting "motion" of the carriers in the energy space and it is valid when r(p) << t « te, 
where r(p) is the momentum relaxation time and te is the energy relaxation time. The method is semiclassical by 
its own nature and applicable when the energy exchange between the carriers and the surrounding medium can be 
assumed quasicontinuous which excludes highly inelastic scattering processess. It is valid when the average carrier 
energy is much larger than the exchanged energy as in the case of high-field transport. The FPA has the advantage of 
being analytical, and, whenever it can be applied, saves computational time and allows a more transparent physical 
\ interpretation. 

In this report, we present results within the FPA for bulk Si and compare them with experimental data and previous 
Monte Carlo simulations (using the BTE). We show that the FPA leads to good results as compared with those data 
whatever several scattering mechanisms (six intervalley phonons between the Si A valleys in the conduction band) 
are taken into account. The intravalley acoustic phonons were ignored, and thus our results are reliable just for high 
enough temperatures. 

The evolution of the energy distribution function (DF) f{E,t) is governed by the Fokker-Planck equations 



where 



J(E, t) = W(E)N(E)f(E, t) - A [D(E)N{E)f{E, t)} , (2) 

such that f(E,t)N(E) gives the number of carriers at time t with energies in the interval (E,E + dE), while the 
function N(E) represents the density of states (DOS). In Eq. (||) W(E) represents a certain "drift velocity" in energy 
space and in fact gives the rate of energy balance of the carrier, D(E) is a kind of diffusion coefficient and J(E,t) 
represents thus the carrier current density in energy space. Under steady state conditions, the Eq. @ transforms 
into 

A [D{E)N{E)f{E)\ = W(E)N(E)f(E). (3) 

The carriers interact with the phonons and the applied dc electric field F. We assume a phonon reservoir in thermal 
equilibrium at the temperature T and that the continuous exchange of phonons between the carriers and the bath 
does not affect the thermal equilibrium of the latter. Hence, the coefficients W(E) and D(E) are split as follows 

D(E)=D F (E) + D ph (E) , W(E) = W F {E)+W ph (E). (4) 



1 



The label "F" ("ph") denotes the electric field (phonon) contribution to these coefficients, whose explicit forms will 
be given below. Equation (^|) has the simple solution 

f(E) = cx P jy 

where D p h{E) was neglected. This approximation is very well fulfilled in all the cases of interest for us ||. 

We consider transport of electrons in the Si conduction band (CB) in a high dc electric field regime (F > 10 kV/cm). 
We take into account the six ellipsoidal energy valleys of Si at the A points of the Brillouin zone (along the < 100 > 
direction). To be specific, let us take F = (0, 0, F), where the z axis is taken along one high symmetry direction, and 
denote by "Z" ("tr") the valleys with principal axis parallel (perpendicular) to F. 

The explicit expressions for Dp(E) and W p h(E) are 

D F (E) = (r(p}[ q F-W f e(p)} 2 ), (6) 



WphjE) 
D F (E) 



dE 



(5) 



W ph (E) = %UJ [l/T abs {p) - 1/Tem(p)] , (7) 

where the brackets represent an average over the constant energy surface e(p) = const, r(p) is the total relaxation 
time due to the electron-phonon scattering, and "a&s" ("em") denotes phonon absorption (emission) by the electron 
due to the several scattering mechanisms. 

A straightforward evaluation of Eq. (^J) leads to 

D F(E) = ^^r(e) 1 (e)/( 1 '(e)) 2 j = l,tr, (8) 

where j'(e) denotes the derivative of the function 7(e) = e(l + ae) responsible by the non-parabolicity of the band 
structure with e = e(p) being the energy dispersion for a given valley. For r(e) we shall consider the six intervalley 
phonons responsible for transitions between the equivalent A-valleys of Si. Then the total relaxation time reads as 

1 6 r 
-^^VOn m(T) ^7(e + + 2a(e + h^)\ 

r(e) ^ L 



+{ni{T) + l)V7(e-^)|l + 2a(e - hui)\Q(e - , (9) 

where O(e) is the step function, rii(T) is the phonon distribution function, and C i = (m tr m]^ 2 D'^ i )/('\/2T:ph 3 LUi) with 
Wtr (™z) being the transverse (longitudinal) effective mass, p the semiconductor density, u>i and D i are the frequency 
and deformation-potential constant respectively for intervalley phonons of type i || . In Eq. (^|) , the first and second 
terms correspond to 1/ T a& s ( e ) & n d 1/Tg m (e) respectively. The phonon contribution can be written as 
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Wph (e) =Y / W; h (e). (10) 

i=l 

where W^ h {e) is given by Eq.(Q) for each i. Intravalley optical phonons do not contribute to transition rates, because the 
corresponding transitions are forbidden by the selection rules and we assume that, for high T and E, the contribution 
of intravalley acoustic phonons should be neglected. 

Once, we have evaluated W p h{E) and D F (E), we have to calculate the integral in Eq.(|5|) to obtain the two DF, 
fi(E) and ft r (E), corresponding to the l-tr valleys respectively. Of course, the FPA is of practical use only in the case 
when such integration can be analytically performed which is not the case for expressions of D F and 1/t(E) given by 
Eqs.(|) and (§). So, we consider a single effective intervalley phonon with energy Tilhq = 0.0343 eV, obtained from an 
average of different phonon frequencies given in Table VI of and with constant D a obtained by the superposition 
of the different deformation-potential constants, but also including the number of final equivalent valleys for each kind 
of transition. With this approximation we obtain the following result 



DUE) ~ 2e 2 F2 



Y J Co l C ok <$>{E,T), j = l,tr, (11) 
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with 



$(E,T) = n 2 {T){E + 1)(1 + E {E + 1))(1 + 2E (E + l)) 2 

- (n(T) + 1) 2 (£ - 1)(1 + £ (£ - 1))(1 + 2E (E - l)) 2 G(E - 1). (12) 

where hereafter E is in units of Huq and Eq = HujQa. Considering the contributions from different valleys and using 
Eq.(||), again applying the parameters from Table VI Ref. ||, we are led to: 



f(E) = exp 



(3j / <£(x,T)dx 



(13) 



with f3 j = (?>mjrn d D A )/ATT 2 p 2 e 2 F 2 fi A and D Q is the effective deformation-potential constant defined through D A = 
D 2 jD 2 k . We have estimated D — 12.09 x 10 8 eV/cm and in all the above expressions the overlapping integral (see 
Ref. Pj) was taken equal to unity. However, for numerical computations it should be useful consider it as a fitting 
parameter. 

The integral involved in Eq. (|l^) can be analytically performed in a straightforward way and the general structure 
of the DF has the form: 

f(E) = E A (l + E E) B exp(C-P(E)), (14) 

where A, B and C are parameters dependent on T and F and P(E) is a polynomial.. This structure is far from a 
Maxwellian one. The DF describes the stationary non-equilibrium configuration where an electron temperature T e 
cannot be defined. From Eq.(jl4|), we can immediately obtain the average electron energy, estimated as E av (T,F) = 
(E l av + 2E%)/3 with 



El 



av 



Efj(E)N(E)dE 



I \ f(E)jN(E)dE j = I, tr. (15) 



In Fig. 1, the electric-field dependence of E av is depicted for different temperatures. As expected, we found a weak 
temperature dependence. We can see that the average electron energy increases for larger electric fields. Moreover, 
the condition E av >> Tiujq is fairly well accomplished, ensuring that the FPA is within its range of validity for the 
given temperatures. Our results cannot be expected to be correct for low temperatures (or too low carrier energies) 
because we neglect intravalley acoustic phonons. 

The drift velocity va = (v<u + 2v,itr)/3 can be also evaluated from |||| 



_ 2eF f 1 (E)t(E) 
Vd3 ~3^J (HE)) 2 



dfj(E) 



dE 



N(E)dE/ / f(E)jN(E)dE. (16) 



In Fig. 2, we show as a function of F for two different temperatures. We see that the general behavior of 
the curve is qualitatively correct if compared with experimental results and Monte Carlo simulations. For a more 
quantitative comparison, we present in Fig. 3 our results together with the experimental data and those from Monte 
Carlo simulations |l^Jll[]. As it can be seen, we obtained a good agreement with both results for T = 300 K in 
the high electric-field regime. However, in the opposite limit, the results from FPA do not reproduce those from 
experiments and Monte Carlo calculations, as should be expected. For comparison with experimental data we have 
taken j3j as a fitting parameter, an issue which can be reasonably understood considering the overlapping integral for 
the electron-phonon scattering probabilities. Another point to be stressed is that reliable results were obtained just 
for high temperatures. However this is not a limitation of the FPA itself but of our present calculations since we have 
ignored intravalley acoustic phonons. 

In conclusion, we have shown that transport problems can be tackled by the mathematically simple FPA even in 
the case of a concrete semiconductor. We can also emphasize that the possibility for achieving correct results from 
the FPA depends on the chance of performing good enough approximations for the relaxation processes, allowing the 
analytical evaluation of the integral in Eq.(^|) and still retaining the essential physical picture. The results, being 
acceptable just for the higher electric fields, in fact are very close to those of Monte-Carlo simulations. Additional 
calculations for T= 430 K also revealed good agreement with those of [0. The saturation effect, however, is not 
predicted by FPA. For higher electric field intensities a slow but ever decreasing behaviour is achieved. 
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FIG. 1. Average electron energy (in units of Tiujq) as a function of the electric field for three different temperatures. It is 
clearly seen that the average energy is actually much larger than the phonon energy. 

FIG. 2. Drift velocity for three different temperatures as a function of the electric field. The general trend of the curves 
corresponds with what is seen in experiments and Monte Carlo simulations. 

FIG. 3. Drift velocity as a function of electric field for 300 K for Si. Our results (solid curve) are compared with Monte 
Carlo simulations (dotted curve) and experimental data . 
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